ivrfvio.ml Advisory Committee, 
tor Aeronautics 
MA ILED 

APR 171939 

i Library, 



TECHNICAL NOTES 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



"TVl- ~ft>0 

No. 700 



THEORY OF. AUTOMATIC CONTROL OF AIRPLANES 
By Herbert K . WeiSB 



\j t . .-. a ■. ■ RESEARCH GEfcttES 
-sARY.NASa "i^r 
(.tn^lON t VIRG!llia* 



la. shin gt on 
April 1939 



' • A 



v 9 



NATIONAL ADVISORY COMMITTEE EOR AERONAUTICS 



TECHNICAL NOTE NO. 700 



THEORY OE AUTOMATIC CONTROL OP AIRPLANES 
By Herbert K , Weiss 



SUMMARY 



Methods of, automatically controlling the airplane are 
reviewed. Equations for the controlled motion including 
inertia effects of the control are developed and methods 
of investigating the stability of the resulting fifth and 
higher order equations are presented. The equations for 
longitudinal and lateral motion with "both ideal and non- 
ideal controls are developed in dimensionle ss form in 
terms of control parameters "based on simple dynamic tests 
of the isolated control unit. " 



INTRODUCTION 

Automatic control implies the process of making some 
physical quantity take on an arbitrary and predetermined 
series of values without human supervision. A perfect 
control for aircraft would maintain the airplane along a 
desired flight path and would completely suppress undesired 
disturbances in pitching, rolling, and yawing. 

The means for applying such a complete constraint to 
the airplane are lacking.. In the conventional airplane, 
the pilot can influence the motion only "by movement of the 
elevators, the rudder, the ailerons, and the throttle. 
The law by which these controls are adjusted can be related 
to any characteristic of the motion, but the controlling 
influences can be applied only as rolling, yawing, and 
pitching moments, and as a longitudinal force. 

The problem of automatic control lies in relating 
these controlling influences to the natural characteris- 
tics of the airplane so as most nearly to attain the per- 
formance of the perfect control. 

In order to study the motion of the controlled air- 
plane, it is necessary to extend the equations of motion 

*Basedjan thesis submitted in partial fulfillment of the 
requirements for the degree of Master of Science, 
Massachusetts Institute of Technology. 
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to include the physical characteristics of the control, 
which in any actual case will not act instantaneously and 
may possess various kinds of lag. In the following paper, 
therefore., atte-ntion is first given to the performance of 
a generalized control, isolated from its controlled mem- 
ber, and the complete controlled motion of tho airplane is 
then established in terms of parameters based on the free 
motion of the control. 



STUDY OP THE ISOLATED CONTROL, WITH METHODS FOR 
DETERMINING THE STABILITY OE HIGHER ORDER EQUATIONS 
Control Charaot oris tics 



The controls to be discussed are all "error sensi- 
tive"; that is, they operate to maintain some quantity 
constant but derive the impulse for their operation from 
an error in this quantity. Wiiile they can make the error 
very small, they cannot entirely eliminate it. 

Three degrees of sensitivity to the error may bo notod. 
The control force may be a function of the error magnitude, 
of the rate at which the error is changing, or of the sec- 
ond derivative of the error. In the most general case, the 
controlling force is proportional to both the error and its 
derivatives, and may be expressed as 

E=a 1 <+a s j+ a 3 e (l) 

Minorsky (reference l) has suggested that the error 
and its derivatives might also govern the rate at which 
tho controlling force was applied. The two additional 
cases that he advances may be written as 

dE/dt = bic + b a c + b 3 e (2) 

I ... 

d 2 E/d 2 t = c 1 € + c 2 e + c 3 e (3) 

whore 

£ is tho error (differenco botwoen tho valuo 
desired and the actual valuo of tho con- 
trolled quantity). 
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F, the controlling force.- ■' - 

a, b, and c, constants of proportionality. 

The three types will he called Class I, Class II, and 
Class III, respectively. It will he seen that Class I 
controls -allow a constant error when there is a steady 
disturbing force. This error can he reduced hy increasing 
the control sensitivity, "but it cannot he completely elim- 
inated. Most automatic controls for airplanes fall in 
Class I and are usually without "benefit of the derivative 
components (a a and a 3 =-0). 

Case II controls are used when the controlled quanti- 
ty is subjected to prolonged and slowly changing disturb- 
ances. They allow no steady-state error under constant 
disturbing influences for as long as an error is present, 
the controlling force increases. Glass III controls admit 
steady-state errors only when the second derivative of the 
distur hance varies. 



The Control as a Simple System 

In most cases, the assumption can he made that the 
control is equivalent to a simple system with only One de- 
gree of freedom. The exceptions are the controls in which 
two of the components have approximately equal frequencies; 
in this case, the isolated control may develop peculiari- 
ties corresponding to a system with more than one degree 
of freedom. ^ for practical computations, however, the con- 
trol may "be replaced hy an equivalent mass., equivalent 
inortia, equivalent "static stability, " and coupling ratios 

The control inherently possesses inertia. The smaller 
its inertia, the more satisfactory a control is likely to 
be because the -"inertia lag" is reduced. Damping is often 
added by design to eliminate the tendency of the control- 
controlled system to hunt. Servo mechanisms embodying hy- 
draulically operated pistons may possess the equivalent of 
damping because of the resistance of tho fluid in the sup- 
ply lines to change of velocity of flow. 

"Static stability" of the control .require s that a 
small departure from the neutral position should produce 
a. force in the control tending ,to return the control to 
the neutral position. When the -e^ntrollod quantity has no 
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inherent stability of its own, the control trill generally 
require static stability. Thus, an azimuth control for 
airplanes must he statically stable, in order that . the 
airplane -» itself insensitive to direction in azimuth - 
may hold a given course. 

It will he assumed that the controls discussed are of 
the type which Hazen (reference 2) calls "continuous con- 
trols." Any error, however small, is considered to pro- 
duce a corresponding controlling influence -through the 
mechanism. Actual controls may have a small inactivo zono 
within which they are insensitive to errors. The motion 
of these controls can he determined by solving for the in- 
active and active regions separately, with due considera- 
tion for boundary conditions. 



Response of the Isolated Control 

The dynamic characteristics of a control isolated 
from its controlled system can be obtained by subjecting 
the control to an arbitrary forcing function. The sim- 
plest disturbance consists of the sudden application of a 
constant force or displacement to the control. The theory 
for these "step functions" will now bo devolopod. 

Lot 

x be the position of some characteristic point on 
the control, referred to the neutral position 
of the control. 

m f equivalent inertia of the control referred to x. 

c, equivalent viscous resistance (damping) of the 
control referred to X . 

k , static stability of the control referred to x. 

L, a load suddenly applied to the control. 

x^ , a displacement suddenly applied to the control. 

Coordinates may, of course, be linear or polar, de- 
pending on the physical arrangement of the control. 

The equation of motion is, when the force step func- 
tion is considered, 
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L - mx - cx - kx = 0 ' (4) 

' Then let 

uu n "be natural angular . frequency of the control ( </ k/m ) . 

£ = damping ratio, the ratio of the amount of damping 
present in the control to the amount necessary to 
produce motion that is just short of being oscilla- 
tory.. In terms of t and uu n , equation (4) is 
then 

X + 2£<%x + to n 2 3C = J u^ 3 (5) 

and the steady-state displacement of the control is 

*.« - L /k (6) 

The general solution for oscillatory motion (£ < 0) 
may "be developed as 

x/x gg = 1-e^* (co s yrrpou n t + — sin^/l - t a u) n t) 

v 1 — i 

If the motion is undamped, this equation simplifios to 

x / z ss = 1 ~ cos lUa* ( 8 ) 

and, when £ = 1.00, the motion is critically damped and 
is expressed by 



:/x ag = 1 - «-«%* - u> n t-^t " ; " ( 9 ) 



A. damping ratio greater than 1.00 produces motion 
similar -co critically damped motion but more sluggish and, 
as it is essential that the control should operate quickly, 
the overdamped case will not be of great importance. 

Figure 1 illustrates the effect of varying the amount 
of damping on the response of the control. 

When the control is suddenly displaced rather than 

being disturbed by a force, the o quaTl ins express sin ? its 

return to equilibrium are the sane aB those presented with 
the exception that tho initial 1.00 is lacking. Thus, if 
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an aperiodic control (equation (9)) Is displaced so that 
its position of equilibrium varies "by x^. , the control 

attains its new equilibrium according to the expression 

= -c^te^ (10) 

x^ ■"• 

due attention being paid to the proper signs. 

Rapidity of response is desired and, since the motion 
becomes slower with increased damping, it is evident that 
from this consideration damping should be small. Suffi- 
cient damping must be retained, however, to insure that 
the corresponding oscillation rapidly decays. 



Lag in the Control 

For most controls it is a straightforward process to 
set up the equations of controlled motion, assuming that 
the control has no inertia, dead period, or friction, and 
acts instantaneously. 

It is customary to lump the deficiencies of actual 
controls that prevent them from attaining this "ideal" per- 
formance under the general heading of lag, and to calcu- 
late their effects approximately by some semiempir ical 
method. 

Strictly speaking, a distinction should be made be- 
tween the lag of an inactive zone at the neutral position 
of the control, and the lag that extends over the whole 
range of operation of the mechanism, such as that caused 
by the inertia of the parts. 

The methods that have been used in treating lag in a 
control may be grouped roughly under f-our headings: 

1. Introduction of inactive zone of control. 

2. Assumption of constant time lag. 

3. Use of semiempirical approximations. 

4. Use of control characteristics. 

Method 1 has been demon-st rated by Eazen (reference 2) 
and Klemin (reference 3). The solution for the motion is 
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obtained in parts, with, due consideration for the boundary 
conditions "between regions of active and inactive control, 

T»hen lag is known to "be present hut cannot he exactly 
evaluated mathematically, as in the case of .lag of the hu- 
man pilot, several writers have used the approximation 
that the control response lags the error hy a constant time 
interval. Minorsky (reference l), Oallendar, Hartree and 
Porter (reference 4), and Cowley (reference 5) trerit lag 
hy method 2. 

Minorsky uses a Taylor's series directly so that, if 
the controlling influence is a function of conditions m 
seconds previous, it may he expressed in terms of present 
time t as 

f(t-m) = f(t) - mf'.(t) + ffl 8 /21 f"(t) - ... (11) 

and if m is small, higher order terms can he neglected. 

Oallendar, Hartree, and Porter approach the prohlem 
hy making the usual assumption that the solution of the 
differential equation is of the form 

x t = Ae** (12) 

and then, if the control moves according to x at time 
t-m., the terms expressing the control influence depend on 

x^ = Ae X(t ~ m > (13). 
t— m 

The resulting equation is no longer linear in 7\ 
hut can he solved hy expanding the exponentials in a se- 
ries and neglecting higher order terms, or "by graphical 
.means (reference 4). 

Garner (reference 6) has used a simplified method of 
treating control lag with empirical constants, which 
amounts to a consideration of the first two terms of ei- 
ther of th.e foregoing series. The semiempir ical approxima- 
tions of method 3, while useful for a rough check of the 
effects of control lag, have the difficulty of employing 
arhitrary constants not always availahle in any particular 
case. 

Method 4, the introduction of the over-all frequency 
and the effective damping of the control into the expression 
for the controlled motion, has "been followed in the present 
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paper. This procedure has the advantage of using quanti- 
ties that can "be experimentally measured in the laboratory. 

The introduction of additional - degree s of freedom into 
the already, complex expression 'for the motion of the air- 
plane necessitates the development of-a method for the 
treatment of fif-th and higher order equations . This meth- 
od will now be considered. 



Higher Order Equations 

In order that the motion of the airplane lend itsolf 
to mathematical treatment, it is necessary either that all 
relations be linear or "that only small motions be consid- 
ered. The classical treatment of the motion of the air- 
plane, has by now been just if ied 'as applicable to disturb- 
ances of appreciable magnitude. The linearity of control 
response depends on the design "of the control, but violent 
movements will not be expected and a linear response for 
small displacements is a fair assumption. ... 

Solution is obtained by' writing the differential equa- 
tions of motion, assuming a solution of the form 

x = S A k e*k* . (14) 

and expanding the resulting determinant into an equation 
linear in A. The difficulty lies in the solution for the 
various values of A from the equation, which is of the 
form 

aA a + bA 11 - 1 + cX n - s + dA n ~ 3 +.•..= 0 (15) 

Methods are available (references V to 13) for the so- 
lution of the quartic equation and for expressing the com- 
plete motion of the uncontrolled airplane. There are also 
methods of extracting the xoots, complex or real, of the 
quintic, sextic, and higher order "equations (references 14, 
15) that result from controlled motion; but the methods 
are long and troublesome, especially when it is desired to 
investigate a range of possible variations of the control 
relat ions . 

It will often be sufficient to determine ■ simply tho 
stability of the controlled "motion. * 
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Methods of South and Hurwitz 

South (reference 16) presents a series of tost func- 
tions which can "be built up for any decree of linear equa- 
tion and which indicate from the coefficients of the equa- 
tion whether the motion it represents will increase or de- 
cay with time. The fun'ctions are obtained by. writing the 
sequences 



a b c d'\ <a~, . £\ SJ. 

h ("bc--ad) d (be-af) f (bg-ah) h 



(16) 
(1.7) 



Beginning with a, each. test function is derived 
from the one preceding it by substituting for each letter 
in sequence (16) , the letter or expression directly below 
it in sequence (17.) (substituting zero for letters above 
those appearing in the original equation)^ The motion" is" 
then stable if the final function and the coefficients are 
positivoi 

Obtained in this manner, the Routhian test function 
for the quintic is 

I (be - ad) d - b (be - af)l (be - af ) - (be - ad) ^f 

J (18) 

Hurwitz (reference 17) giye.s a method of obtaining 

the stability functions as an expansion of the determi- 
nant 
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The-motion is stable if the determinant and the coef- 
ficients are positive. .This f.orjn is somewhat simpler than 
the one used by South for numerical substitution but still 
involves considerable work far equations above, the quartic. 

Sefore a simpler method of following, stability changes 
as the relations of a mechanical system are systematically 
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varied is presented, the form of the solution of the oqua- 
tions of motion will be discussed in more detail. 



Form of the Solution 

When some of- the values of A-^ to be substituted in 
equation (14) are complex (and therefore conjugate pairs 
since the coefficients of the original equation were real), 
there are pairs of roots of the form a ±ib. The part of 
the solution corresponding to these two roots can then bo 
written 

x = e at (A n cos bt + A ai sin bt ) (20) 

The test functions of Houth and Hurwitz, when posi- 
tive, insure that a be negative, so that the oscillation 
represented, decays with time. If a is zero", so that K 
is a pure imaginary, the term represents an unending os- 
cillation, which is the boundary condition between stabil- 
ity and instability for the term. 

The presence of pure imaginary terms is indicated by 
the fact that South' s discriminant becomes zero. It should 
be noted, howeve-r, that the discriminant is also zero for 
more than one pair of equal pure imaginary roots and, when 
two sets of roots are ^i.a = ^3.4 = ± ib, the motion is 

unstable, being of the form 

x = Ax cos bt + A s sin bt + A 3 t cos bt + A 4 t sin bt 

(21-) 

Transition to Instability 

The fact that a pair of roots becomes pure imaginary 
as the system from which the equation is derived passes 
from a stable to an unstable condition ib made use of in 
determining the point of critical stability when some 
physical characteristic is varied systematically. 

Define the angul ar frequency of a complete mechan- 

ical system of any number of degrees of freedom as the fre- 
quency at which it can execute unending oscillation. Sim- 
ple systems of one degree of freedom will oscillate endless- 
ly only in the absence of damping, but more complicated 
systems easily ana sometimes annoyingly perform self-excited 
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oscillations in spite of a large amount of damping (refer- 
ences 18, 19). .The expressions for natural frequency will 
therefore he. estahlished in their most general form, in- 
cluding damping. The most general case of unending oscil- 
lation for the systems that can support more than one end- 
less motion will he considered as the case in which the 
corresponding equation has only one pair of imaginary roots. 



Expressions for Natural Frequency 

At the "boundary of transition from stability to in- 
stability, it is known, that two of the roots of the equa- 
tion in A are ±iuu n . Either root can then he substi- 
tuted "back into the original equation, the sum of the real 
and the imaginary terms he equated to zero, and tho reduced 
equations solved for tho frequency. Performed in detail 
for the cubic equation, the process is as follows: 

ax 3 + bx a + ci + d =. 0 

-iauj n 3 - bu) n a + cit% + d = 0 

tu n s , * c/a t (22.1) 

oi n a o d/b • (22.2) 

Unending oscillation is indicated when numerical sub- 
stitution produces the same value of U) n 3 from both ex- 
pressions. It is evident that this method is equivalent 
to setting Routh r s discriminant for" the cubic (be - ad) 
equal to zero and, for this simple case, there is no gain 
in simplicity by performing the operation in two parts. 



Eor the higher order equations, the pairs of expres- 
sions in tu n are, when similarly obtained, 



QUART I C 


U)„ a 
n 


= d/b 






(23.1) 




< 


= ( cd - 


be )/ba 




(23.2) 


QUINT IC 




= (be - 


af )/(bc 


- ad) 


(24.1) 






= (de - 


cf )/(bc 


- ad) 


(24.2) 
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SEXTIC ujl 



2 _ ("bo - ad)f - 13 g 



(25.1) 



("be - ad)d - b(be - af ) 



^n* - (cL/Td ) 0J n 3 ■+ (f/b) = 0 



(25.2) 



SBPTIO (ju. 



a 



("be - ad)(fg n. eh.) - (be - &h ) 



'n 



("be - ad)(dg - ch) - (bg - ah)(be - af ) 



(26.1) 



("be - ad) 0) n 4 - (be - af) cu n a + (bg - ah) = 0 



(26.2 ) 



The expressions were brought into tha form given by- 
straightforward mathematical "Juggling." When a single 
computation of the stability - of a mechanical system is de- 
sired, Hurwitz 's determinant should be numerically expand- 
ed, as this method will be found simpler than the several 
computations necessary to determine, by moans of the nat- 
ural frequencies, whether the system is on the stable or 
unstable side of the critical point. 

J7hen it is desired to observe the transition from sta- 
bility to instability of a mechanical system, as when stud- 
ying the effects of control lag ox inertia, the frequency 
method is much more convenient than that of either Kouth 
or Hurwitz. In addition, when the transition point has 
been determined by this method, the frequency of the end- 
loss oscillation is at once available, without further 
computation. 



An example of- the application of tho method to a 
practical caso will now be borrowed f_roin a later section 
of this report. It is desired to learn within what range 
of free naturals period a longitudinal control of a givon 
sensitivity will bo satisfactory, i.e., will not allow 
s-olf-excited oscillation. 

The equation involved is a sexiic. It has been de- 
termined by applying Routh's discriminant to the quartic 
equation for the airplane motion under ideal control that, 
the motion is stable when the natural period of the con- 
trol is infinitely short (T^ - 0). As soon as a finite 

natural period is admitted, the equation of motion becomes 
a sextic. In figure 2 the two expressions for natural 



Application to Sextic Squation 
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frequency and Routh ' s discriminant .for the sextic have 
"been calculated and plotted against increasing natural 
period of the control. 

By either method, the motion is seen to "become unsta- 
ble when the free natural period of the control exceeds 
1.95 seconds, hut the actual computations were much sim- 
pler for the frequency curves. The frequency of the un- 
ending oscillation is, "by inspection, ahout 19 radians per 
second. It is interesting to note that the frequenc y of 
the short and heavily damped oscillation of the ideally 
controlled airplane considered is 16.9 radians per second. 

When the system is nearly critical, it is possible to 
calculate the rate of growth or decay of the nearly end- 
less oscillation by a method of Blondel's (reference 20) 
which is based on the assumption that a in the root 
a-ib is so small that higher powers of a are negligible 
compared with a itself. The root is then substituted 
back into the original equation in >», and real and imag- 
inary terms are separately equated to zero, when solution 
can be made for a and b. 



AUTOMATIC CONTROL 0* AIRPLANES 
LONG- ITTJD INAL MOTION 



If the airplane is slightly disturbed in smooth air 
and allowed to execute free longitudinal motion, it will, 
if dynamically stable, regain a steady-flight condition as 
the disturbed motion decays in the form of two damped os- 
cillations. These two modes of oscillation consist of: 

1. A heavily damped oscillation of short period (of 
the order of a few seconds) involving primarily change of 
incidence, in which changes of forward velocity are negli— 
gible-.- This motion disappears almost at once, and in most 
airplanes is not noticeable as an oscillation. 

2. A long period, lightly damped oscillation involv- 
ing change of forward speed, during which the airplane 
rises and falls. This oscillation depends on the drag of 
the airplane for its damping, and is increasingly trouble- 
some on "clean" airplanes. 

Although the short oscillation has been somewhat neg- 
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looted in stability analysis, "because it disappears so rap- 
idly, Jones (reference 21) po int s . out . that , although the 
heavy damping of this mode insures its rapid subsidence in 
calm air, it imposes an effective restraint against move- 
ments of the airplane relative to the air and causoB vio- 
lent movements of the airplane in- rough air. 

Law of Operating Control 

The controlling moment in longitudinal control — is ex- 
ert-ed "by means of the elevators, which are then to be moved 
according to some flight characteristic. Haus (reference 
22) >ives the foil owing table of disturbance detectors and 
the quantities to which each is sensitive. 



In strument 

1. Air-speed indicator 

2. Wind vane 

3 . Free gyroscope, sus- 

pended at its e.g. 

4. Motor-driven gyro- 

scope with proces- 
sional moment 

5. Pendulum or acceler- 

ometer along OS 

6. A»c c elerometer along 

OZ 

7. Lift indicator 

8. Rate-of-climb meter 



9. Tor sional- accelerom- 
cter about OY 



R ecorded quantity. 

Relative speed 

Incidence 

Absolute inclina- 
tion 

Angular velocity 



Direction of appar- 
ent gravity 

Magnitude . of- appar- 
ent gravity 

Magnitude of lift 

Speed along verti- 
cal 

Angular accelera- 
tion 



Sym bol 
U 

a = -w/u 
6 



du/dt and 
r, in 6 

dw/dt and 
. cos 8 

iV s or uw 

w ,or V sin 6 

6 



The elevators can be moved according to the indica- 
tions of ,-,ny of those instruments or combinations of them. 
The most successful controls, those of Sporry (references 
23, 24, 25) and Smith (references 26, 27) are of type 3, 
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operating according to the absolute inclination -of tho air- 
plane in spaco . This type '•will 'therefore "be considered in 
detail. Analysis of tho other control types would ho car- 
ried out in a manner similar r to the one now presented. 



Equations of Controlled Motion 

In its broadest form, the type. 3 control possesses ve- 
locity and acceleration -components (types 4 and 9). A gen- 
eral solution including the displacement-, the rate, and the 
acceleration components is .no more difficult than that for 
the simple control, and .the .full form will thereforo be 
considered. 

The elevators are linked to tho control through a 
servo mechanism, so that the pi t.ching moment varies accord- 
ing to the displacement, of -the control. For small motions, 
a linear relation can bo assumed and, for many controls, 
the assumption will also be valid for large displacomont s . 

Since the control mechanism has inertia, an additional 
degree of freedom is introduced, and there are now four si- 
multaneous equations of motion: 



X 


= m(u 


Wq.) 


(27) 


z 


= m(w 


- Uq.) 


(28), 


M 


= b4 




(29) 

■ . ....... ^ 


F 






• (30) 



where P is force on the control. 

m , effective inertia of the -control referred 

C tO £ . : 

|, displacement of the control. 

The other symbols have their usual significance. 

The full expressions for horizontal and vortical force 
are as usual, but equations (29) and (30 ) are now. 



Bq = u M + w M + q M - t M (31) 
. 8u 3w * 8q. b 3| v ' ' 
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m c £ = 6 §5 + * §T.+ 6 fs - 5 ^7 - t ^~ 



And, of course, % 



m 



a5 ~ b at °.a| 

is the controlling noment 



(32) 



corresponding to a control displacement |, and 



b 33 



is the force tending to operate the control for an angular 
displacement of the airplane 6, The remainder of the de- 
rivatives signify corresponding linkages. 

Writing now the full determinant for the complementary 
solution where D is the operator d/dt , 



D - X, 
0 



u 



w 



-g 



33 Z T 
■ -1L 
0 



-gV - DU 0 
D 3 - DM a 
-D 3 Fg -DP -Pg 



0 
0 



(33) 



It has "been assumed that most of the inertia of the 
control is effective after the three impulses have been 
combined to operate the elevators, which is equivalent to 
assuming that the error and the derivative controls have 
equal effective inertias and dampings. While experimental 
tests are required to determine the accuracy of this assump- 
tion, it is probably adequate for well-built controls, and 
controls could certainly be built for which it would hold 
exactly. 

Note that in equation (33) the minor consisting of. 
the first three rows and columns is the determinant for the 
uncontrolled airplane. Set the whole determinant equal to 
zero for the complementary solution and expand it in torma 
of this principal third order minor. This procedure givos 



D + DP? + 



D - X 



u 



-Z 



u 



-X w 
D - Z 



w 



= 0 



(34) 



This form can be solved, but the final result will be 
more advantageous if it is put into nondimensional form'. 
Going back to expression (33), write the dimensions of 
each term. 
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In terms of its dimensions , the determinant "becomes 



GT 1 


5T X 


LT~ a 


0 


3,-1 


5,-1 


LT""' S 


0 




L~ 1 T - 1 


T -2 


L -i gl-B 


0 


0 


LT~ S 


T~ 2 



Following the same procedure as in nondimensional iz- 
ing the uncontrolled motion of the airplane, multiply the 
derivatives of 

the first row by T 
the second row by T 
the third row "by LT 
the fourth row "by T 
the third column by L X T 
the fourth column "by T 

The characteristic length of the dimensionle s s sys- 
tem is taken as L, the length of the tail moment arm. 
The characteristic time is defined "by 

T = mXp/2 SU) (36) 

On this oasis, the unit of velocity is 

L/T = U/'u. (37) 

where p, = m/(p/2 SL) (38) 

and \x may he called the "relative density" of the air- 
plane, "being proportional to the mass and invorsely pro- 
portional to the cube of the linear dimensions of the air- 
plane. .Glauort notes (reference 28) that apart from the 
derivative coefficients of the airplane, H is the only 
parameter which affects the stability. 

The dimensionless form of the derivatives is as do- 
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fined "by Metcalf (reference 29) and other writers, "but it 
is necessary to. determine the form of the control deriva- 
t ive s . 

In passing, it should he noted. that another dimension- 
less system in use is based on the quantity p instead of 
p/2, as in the present paper. The mathematical expres- 
sions, however, are the same for "both systems except for 
this one difference, and the numerical values of the deriv- 
atives differ only hy a factor of 2. 



Control Derivatives 

It will he shown later that, for the purpose of this 
analysis, it is only necessary to determine the ideal con- 
trol derivatives m Q g. 

W, q. , t? 

The moment exerted hy the tailplane is 

M = C Lt S t p/2 U 2 L (40) 

where 

C Lt is tail lift coefficient. 

, tail area. 

U, steady-flight airplane velocity. 

L, tail moment arm. 

Let 8, elevator angle. 

Following the method of Koppen (reference 30) in non- 
diraen sionali zing the control derivatives, write 

*H - ( dG ^ ) £ s Lu 2 (ad 
se - \rjrJ w 2 s t LU (43J 

where "g"g" * s obtained from wind-tunnel data, and «_ 

depends on the control-coupling ratios and sensitivity. 
Then ' 

U$T S = OM/39) (1/B)T 3 = nm 6 (42) 
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and B = T|mL , where T| is the " mass ■ distribution fac- 

S 



tor" (reference 29) so that 

m8 = Ti VaTV W s 



Similarly 

V (If) (2) * (I) 



(43) 



(44) 



this factor is written m_„ to distinguish it from the 



natural m a of the airplane. Finally 



/9C Li A (7^ (I S t U * L ) 

r 6 = lu 6 = VTas-y VaSV (n m l s ) 



(45) 



Tho f , f , and m» derivatives can ho obtained in 
8 | % 

a similar manner, when it is desirod to evaluate them in 
any particular problem. 



The Dimensionless Determinant 
The dimensionless form of the determinant becomes 



D - x 



u 



— » z 



u 



D-z 



-m. 



u 



w 

0 



-^(CjY+D) 



D 2 —Dm, 



0 
0 
m. 



D 2 +Df|+f| 



and equation (34) becomes 



A 0 + m. 



D S + Df£ + f| 



D - x 



u 



D-z; 



(46) 



= 0 



(47) 



Call the second minor & x • Divide the numerator 
and denominator of its multiplier by f| . Then 



20 



N.A.C.A. Technical Hote No. 700 



_o 
A 



Now 



B S /f| + Df|/f| + 1 



"9 =111 



If 1 5M 



(48) 



e = ae 



so that 



similarly 



F 6 

M 6 = M e ^ 



m ei = m e f| 



(49) 



The suoscript i indicates that ja has not heen extract- 
ed. This extraction could he very easily accomplished hy 
proper evaluation of the other coefficients, such as making 
H m£ tako the place of nig. 

The numerator of equation (48) can now bo written in 
terms of the ideal control derivatives, that is, the val- 
ues which the control derivatives would have if there were 
no lag in the control. That is, 

f m t \ a 

Equation (48) can then he written as follows: 
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as 



Going.- hack 'to equation (32), the control forces sum 

- mj + 6 || + q |I + 6 |f - J -9E - * ^ (52) 

and, when the error impulses are zero, the equation for 
the control alone is 



m 



° !+ *f! + 1 St = 0 (53) 



Compare this form with equation (5) for a generalized 
system of one degree of freedom. 

I + 2 £ w n \ + 0J n a 1=0'. (54) 

so that the F derivatives can be written in terms of tho 
natural frequency and tho damping of 'the isolated control. 
That is, 

F| = a) n 3 • = (2-rr/T n ) s (55) 
E{ = 2£ otfe »-2£ (27t/T b ) (56) 

How f^ = T a ', where T is the characteristic 
time of the airplane, and similarly 

so thon 

f| = (2tt T/T n ) (57) 
f£ = 2t (2tt T/T n ) (58) 

This form is convenient because T n and. £ for tho 

isolated control can bo obtained by simple dynamic tests 
in the laboratory. 

Substituting these values into equation (51 ) gives 

^ 0 -^ 1 (D a m +Dm 2C +p.m e ) + A Q jjD 3 (T n /2TTT ) 2 + 2 t (T^/2-rrT )D J = 0 

(59) 
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This division into three major terms is very conven- 
ient. If the control is nearly ideal, T n is very small 

compared with T, and the quantity in the last bracket 
can be neglected, giving the equation for the ideally con- 
trolled airplane. If the airplane is uncontrolled, the 
second term is zero , leaving the original determinant for 
the uncontrolled motion. 

Although it has not been so noted, during the nondi- 
mon sional izing , the operator D was replaced by th© di- 
mensionless operator DT , which was also written as D. 



Effect of Ideal Control 

Equation (59)is the complete form of the determinant 
for controlled longitudinal motion, an abridged form of 
which, neglecting the derivative components and control 
inertia, is given by Klemin (reference 3). 

The effect of an ideal control on the period and the 
damping of the oscillation will first be dotermined. As- 
suming that £n/T is sufficiently snail so that the third 
term can be neglected, the first two terms of the equation 
can bo expanded into the form, 

aD 4 + bD 3 + oD s + dD + e = 0 (60) 

ffrite f = 1 - mg 

m -i = m + in 
ql q qc 

where ' m^ is the natural damping in pitch of the airplane, 
and m^ c is the effective damping in pitch added by the 
first derivative component of the control. 

It is necessary here to differentiate between the ef- 
fect of m^ and of m^ on the motion. Their effect in 

damping an oscillation once bogun is equivalent. However, 
m^ is derived from the relative motion of the air and the 

tail surface primarily, whereas m- 0 is derived by taking 

mechanically the first derivative of the angle of pitch 
snd using it to operate the elevators. Consequently, an 
airplane with large natural m^ may be expected to exe- 
cute violent motion in rough air, in consequence of the re- 
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straint against motion relative to the gusts. On the other 
hand, an airplane heavily damped in pitch by means of the 
derivative control is restrained not relative to the air 
hut to a set of fixed axes in space, determined by the 
gyro s . 

Protection from gusts, as far as the rotary damping 
is concerned, would then consist of replacement of the nat- 
ural damping by artificial damping relative to space axes, 
as far as possible, with the limitation that the airplane 
must still be controllable manually in the event of fail- 
ure of the automatic pilot. 

Returning to equation (60), the coefficients are 

a = f 

b = -f(x u +z w ) - m ql 

d = -m ql (x u z w -x w z u )+p,m w (x 11 -0 Ii 6)-p,m u (x w +0 L )+p,me (x u +z w ) 

e = ~C L y,m w ( Zu -x u 6) - C^rn^ <x w 6-z w ) - nm e (x u z w -x w z u ) 

(61) 

The effect of m.^ c . upon the damping and the period 

is eqactly the same as the effect of m^ , and so it will 

be sufficient to note here that increasing the sensitivity 
of this component of the control will, over the "normal- 
flight range, increase the damping and lengthen the period. 

It is interesting to note that by making mg = 1.00, 
f can be made zero. This method is equivalent to giving 
the airplane zero inertia in pitch. The quartic equation 
then reduces to a cubic. 

As the mg, ■control is the most widely used, it will 
be considered first. 



Effect of Simple Displacement Control 

Let m*g = m = 0. Lot the coefficients of the quar- 
tic represented by A c be a Q , b Q , c 0 , d 0 , o 0 . Then the 
effect of adding the ra^ derivative is to increase these 

coefficients so that 
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a = a 0 

b * b 0 

c = °o - ^ m G 

d = d 0 + l^m 6 (x u + z w ) 

e = e 0 - M < ra 0 (x uZw - x w z u ) (62) 

The effect of small values of mg on the long and 

the short oscillation has "been calculated for two typical 
airplanes. Airplane 1 is the transport considered by 
Metcalf (reference 29) and Airplane 2 is a small 60-horse- 
power parasol monoplane treated by Soule and Wheatley (ref- 
erence 31). 

Figure 3 gives the variation with m^ of the poriod 

and the damping of the long oscillation. The effect of 
mg on the short oscillation in the range plotted is given 

"by the following table for Airplane 2. 



m Q 


Period 
( sec. ) 


Time to damp to 
^ amplitude 
( sec . ) 


0 


4nr07 


1.08 


-.20 


3.39 


1.19 


-.50 


2.89 


1.40 



The effect of mg on the short oscillation is not 

very great. The period is shortened, as might be expected, 
since mg is a spring constant in pitch. The long oscil- 
lation shortens briefly and then, as its damping increases, 
lengthens its period, and finally "becomes a pair of simple 
subsidences. 

With complete restraint in pitch (mg infinite), the 
equation in D reduces to 



x u z u 



x w z w 



= 0 (63) 
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The success of the simple displacement control in 
practice is explained by the tremendous increase in the 
damping of the long oscillation. The time to damp to J 
amplitude can, in the case considered, "bo reducod from 20 
to 2 seconds. For an ideal control, then, both tho long 
and tho short oscillation would disappear almost at once 
after a disturbance . 

Moreover, the values of mg used in the foregoing 

example have "bean extremely conservative, Klomin (refer- 
ence 3), investigating this typo of control, presents 
Mg = - 2,160 as a practical value for tho airplane he con- 
siders and, in dimonsionloss form, this valuo is oquivalont 
to nig = - 900, approximately. Xlemin finds, in conso- 

quoncc , that the long oscillation disappears completely and 
is replaced hy a pair of subsidences . - 

It should "be understood that the control does not cre- 
ate any new damping. It simply makes a more economical 
use of the damping that is already available in the system. 
The poor distribution of damping "between the two oscilla- 
tions is woll known and the nig control may "bo considered 

as a sort of equalizing valve, which allows some of the 
damping of the short oscillation to flow into the long os- 
cillation. The simple displacement control i3 allowed, 
therefore, "by the peculiar original condition of the sys- 
tem, to produce an effect comparable with that of an error 
and derivative control* Needless to say, this fact is high- 
ly advantageous . from considerations of mechanical simplicity 
of tho control. -• ...... 

Uso of tho Acceleration Component 

The controlling moment can also bo mado to dopond on 
the socond dorivativo of the angular displacement , which 
is oquivalont to increasing or decreasing tho offoctivo 
inortia about tho latoral axis of the airplano. 

The inertia might be increased (m'g negative) to re- 
duce the initial pitching acceleration under the influence 
of gusts, or it might be decreased to permit a more rapid 
damping of the subsequent motion. Simple increase of tho" 
offoctivo inortia will not bo ordinarily toloratod, nowov- 
or, because of its unfavorable effoct on tho stability of 
tho motion. 
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Now eupp'oso a positivo value of m'g woro choson. A 

simple form results if m'g is taken as 1.00, as this 

choice provides that the airplane have zero effective in- 
ertia in pitch. The coefficients of the reduced equation 
then. bear the relation to the original equation - that 

a = 0 



c o "* ^ 3: u z w"* x w z u) = very nearly c 0 

*o 
e = e 0 

A few trial calculations show that m§ affects the 

short oscillation almost exclusively, which might he ex- 
pected on noting that only the first two coefficients are 
greatl,y changed by og. A very good approximate factoriza- 
tion, since b is large., is 

CD + b) (D a + (be - d)/b s D + d/b) = 0 (64) 

and the short oscillation becomes a heavily damped subsid- 
ence , with the long oscillation practically unchanged. But 
tho short oscillation was already satisfactorily well 
damped. Therefore, the use of a positive acceleration com- 
ponent of the control docs not soon to be justified, and 
tho introduction of a negative m'g component would bo sat- 
isfactory only if it woro selective, operating only to op- 
pose movement away from equilibrium. 



Introduction of Control Lag 

Actual controls frequently exhibit a fast residual os- 
cillation, which one writer (reference 27) describes as the 
effect of the control trying to act upon the short oscilla- 
tion. More exactly, this new oscillation is probably 
caused by the additional degree of freedom supplied by tho 
control, and certainly depends on control inertia. 

In order to investigate the nature of this residual 
oscillation, the full form of equation (59) will now be 
used, except that the two derivative drives trill be neg- 
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lected, "because they are not at' present incorporated in 
the Sperry or the Smith controls. The e-q.ua t ion is then 

A D - A x p,m e + A Q j^D 3 (T n /2TrT) 2 + 2£ ( T n /2TrT ) dJ = 0 (65) 
Let 

, r = T /^2-tt . (66) 

Let the coefficients of A q '-"be ' a Q , b Q , c Q , d Q , e Q 

and the coef f ici.ents of A 01 (the ideally controlled de- 
terminant) he a x , b x , c x , d x , e x . 

Then, writing equation (65) as the sextic 

. . aD 6 + bp 5 + cD 4 + dD 3 + jjD 3 4- f D + g - 0 (67) 

the coeff icients, of the- -sextiC: are, ■ in t-erms; of the .coef- 
ficients' of the ideally controlled, and- the uncontrolled 
motion, ». '.".."...' • ' - 



b = 2£r + r 2 b Q 

c = 1 + 2£rb Q + r2.c 0 

d = b 4 +■ 2trc c + r. a d 0 

e = c x + 2-t rd Q + r a e 0 

" f = d x + 2£ re Q 

g = e x - - (68) 

The effect of increasing the natural period and the 
damping of a control for a given v airplane and a given stat- 
ic linkage mV can therefore he carried through in an or- 
derly manner hy -the use of the relations given in equa- 
tions (68). This' procedure has been carried out for a typ- 
ical case. •■ 

Ono of the ideal cases 'treated, was for. mg = - 0.50. 
This case has been extended to include control inertia and 
damping-. "-Tt .is very easy to write- down the sextic, but 
its solution is not an enjoyable task, although there are 
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methods of 1 accomplishing the solution. It is much easier 
simply- to. dot-ermine the effect on the stability of tho mo- 
tion of allowing the control inertia to become finite. 



Effect of Inertia Lag on Stability 
The simple quartic for the uncontrolled motion was 

D 4 + 4.20 D 3 + 11.96 D a + 1.94 D + 1.30- =* 0 

and when mg = -0.50, this expression becomes 

D 4 + 4.20 D 3 + 20.96 D S + 19.40 D + 7.70 =. 0 

This equation is for the airplane of reference 31, 
and the derivatives were based on p instead of p/2. 
The damping ratio of the control was held constant at 
£ ss 0.20 while the natural period T n was increased from 

zero. The stability changes were followed by both Routh's 
discriminant for the sextic, and the use of the simultana- 
ous equations in natural frequency. The results have al- 
ready been given in figure 2. The airplane becomes un- 
stable when tho natural period of tho control oxceeds 1.94 
seconds . 

Effect of Damping- in the Control 

In order to determine the effect of control damping 
on stability, the definitely unstable case of T n = 2.40 

seconds was reconsidered xrlth the damping ratio of the 
control increased to 1.00, that is, critical damping. 
Routh's discriminant then became positive and equal to 
1,275. Increasing the damping then restores stability or, 
in other words, postpones tho critical point of neutral 
stability. 

Although it may seem contrary to common sense to im- 
prove the stability by making the control act loss rapid- 
ly, it will bo remombored that the uncontrolled airplano 
was stable, and adding an infinite amount of damping to 
the control can do no worso than rostoro it to this stato. 
If the airplane were originally unstable, damping might 
eventually havo an adverse effect. 

As in the constant-speed control, where the effect 
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has been discussed in detail (reference 32), damping in 
the control allows the initial surge error to he somewhat 
larger on account of the sluggishness of the control, hut 
successive error surges are reduced because the tendency 
of the control to overshoot has been curbed. 

Damping is not completely beneficial, but it can com- 
pensate for the worst property of controls, which is- iner- 
tia lag. Stability of the motion is attained by damping; 
good performance is attained by reducing control inertia. 



Solution of the Critical Case 

When the natural period of the control. is 1.94 sec- 
onds, the motion becomes critically stable, and an unend- 
ing oscillation is present. From the intersection of the 
two curves for the frequency ,. the frequency of this oscilla- 
tion is uu n s = 19. Two of the roots of tho sextic aro then 

known, and tho sextic can be factored into 

(D 3 + 19) (D 4 + 6.45D 3 + 35D 3 + 37D + 15) =0 

The quartic can be solved by Zimmerman !s method 

(D a + 19) (D s + 1.202) + 0.52 ) (D s + 4.25D + 28.8) = 0 

Compare this with the factorization for the ideal con- 
trolled airplane with no control inertia 

(D a + 1.15D + 0.455) (D 3 + 3.05D + 16.9) =0 

-he long oscillation has hardly been affected. It 
seems Justifiable, then, to obtain a rough factorization 
of the sextic in any case below the. critical by dividing 
through by the quadratic factor of the ideal • quartic , which 
corresponds to the long oscillation. 

The short oscillation is also not greatly modified. 
It is the new, or "residual" oscillation introduced by the 
control that becomes unstable. Because of the characteris- 
tics of the linkages, the residual oscillation apparently 
becomes unstable when its period approaches that of the 
short oscillation. 

The period of the residual oscillation shortens and 
its damping improves as- the inertia of the control is re- 
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duced; the oscillation vanishes completely for an. ideal 
control. 

Determination of Control frequency and Damping 

The convenience of the representation of the control 
as a single degree of' freedom system is seen when the lab- 
oratory procedure necessary to obtain T n , I, and mg is 

outlined. The complete control, with dummy control sur- 
faces suitably weighted to represent their equivalent in- 
ertia, in flight, is mounted on a test platform so that it 
can be easily rotated (in pitch for a longitudinal con- 
trol). The platform is given a very sudden change in in- 
clination, and the response of the control is recorded by 
means of a pencil or equivalent recording means, attached 
to an output member, such as a control push rod of the 
mechanism. Prom the record, which may show a damped oscil- 
lation or an aperiodic approach to the new position, the 
damping ratio' and the natural frequency of the control can 
be determined by elementary vibration formulas. 

The ideal control derivative to.q is determined by 

noting the steady-state control angle for a given angular 
displacement of the disturbance indicator of the control. 
Given this ratio, nig is determined_f rom equation (43). 

Only three quantities are required, and they are the 
three that express the effect of the control on the motion 
of the airplane. Therefore, they also serve as conven- 
ient means of comparing one control against another of the 
same type.' ■ . 

Methods of obtaining advantageous values of the three 
parameters remain in the province of detailed control de- 
sign, but the parameters offer a moans of determining the 
suitability of -an existing control for a particular air- 
plane. 

Suppression of Disturbances in G-usts 

Insuring a rapid decadence, of motions once begun is 
only half the job of a successful automatic control. Tho 
other half . consists of the reduction of initial error 
surges from tho desired course, as the airplane oncountors 
an oxternal disturbance. 

The magnitude of the surge error can be investigated 
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with the aid of operational calculus . or "by means 'of mechan- 
ical methods of solution-, such as the differential analyzer, 

The initial acceleration of the airplane under the in- 
fluence of -various kinds of gust can "be written for each 
of the three degrees of freedom as follows, 



G-ust 



Initial Derivative 
acceleration depends on 



Vertical - w Q x w w o Induced drag 

z w w Q Aspect ratio 

m w w Q Static stability 

Head or Tail - u 0 x u u o Total drag 

z u u o Lift 



m u u Q ' Power application 



Rotary gust - q Q 

z o 

m „q„ Tail size and effi- 

ciency 

*These two are negligible. 

In-order to reduce the effect of the gust in any case, 
the corresponding derivative should "be made small, If the 
derivative is zero, the airplane will not "be affected by a 
gust in that sense. Because the motions in each degree of 
freedom are related, however, an airplane with, say, « 
zero would develop a pitching motion in response to a head 
gust, "but the original forcing function would he appliod 
only as a vertical and a horizontal force. 

The only derivatives that can he modified appreciably 
by the dosigner are the static stability m w and tho damp- 
ing m,,. Experimental, data (reference 33) indicate that, 
in general, airplanes with short period and heavy damping 
do the most pitching in rough air, corresponding roughly to 
large m^. and largo m^. 
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The initial accelerations of the airplane represented 
by m w w 0 and m^<l 0 can then he reduced, and tho stabil- 
ity of the airplane will not ho impaired if- the reduced 
derivatives are augmented hy ta^ and mq C supplied by 
an automatic pilot. The value of mg is not numerically- 
equivalent to m w hut has a similar effect, being, in fact, 
a much more desirable derivative as it allows a more sat- 
isfactory distribution of damping, both m^ and n^ c , be- 
tween the two oscillations. 

Systematic computations by Haus (reference 34) shows 
that, even in the absence of the automatic pilot, reduc- 
tion of m w reduces the violence of the motion. 

Unfortunately, little can be done to reduce the ini- 
tial response of the airplane to a vertical gust. Icr 
the firs-t second or so, th© motion is given very closoly 
by 

w/w Q = 1 - e 2wt/T (69) 

and z w , depending on aspect ..rati*, cannot be reduced ap- 
preciably. This situation might have been foreseen for, 
if the air supporting the airplane rises, the airplane it- 
self must modify its course in space. 

Once started, however, the disturbed motion can bo 
made to disappear much more rapidly with the aid of an au- 
tomatic pilot, as Klemin has shown (reference 3 ) by a num- 
ber of calculations. Klemin also shows that tho vortical 
motion in response to a head gust is very much eased by an 
mg control. Wilson, in an early paper "on the effect of 

gusts (reference .35), also indicated the' beneficial effects 
of a complete constraint in pitch. 

Conclusions 

On the basis of the material presented, somo general 
conclusions can be drawn concerning the operation of a con- 
trol sensitive to tho angle of pitch and its dorivativos. 
In tho absence of complete experimental data, tho conclu- 
sions must bo presented with the reservation that, while 
they satisfy available data, very little data arc available 

1. The control with gyroscopic roforencos givos tho 
airplane sonsitivity with respect to axes fixed in space. 
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2. The simple pitch- control readjusts the propor- 
tions of the system, so that the available damping is more 
equally distributed • "between the long and the short oscil- 
lation, with the result that the long oscillation can "be 
made to disappear as rapidly- as the short oscillation and 
to take on the form of two simple subsidences. 

3. The benefit of the simple partial restraint in 
pitch, therefore, is derived from the ability of the dis- 
placement control to act as an equivalent rato control. 

4. The addition of an actual rate control increases 
the damping in pitch (m^) of the airplane, without in- 
creasing the sensitivity (q. m^) to a rotary gust. 

5. Thore is no advantage in using an unsolectivo 
second derivative control and not much advantage in using 
a selective control. 

6. Inertia in the control introduces a third oscilla- 
tion which can be mistaken for the short oscillation when 
it becomes troublesome. 

V. The third oscillation becomes unstable when the 
inertia of the control is increased beyond a critical value 
determined by the airplane characteristics and the control 
damping ratio. 

8. Damping in the control reduces the effoctiTonoss 
of the control, but stabilizes the residual oscillation. 

'9. Flight in" rough air will be improved (greater 
course stability) by reducing the static stability (m w ) 

and the natural damping in pitch (m^), and by adding an 

automatic control to supply sensitivity to angle of pitch 
(mg) and absoluto damping in pitch (m^ c ). 

CONTROLLED LATERAL MOTION 



The lateral motion of the airplane fo,r small disturb- 
ances consists of translation along the Y axis (sideslip) 
and rotations about the X and Z axes (rolling and yawing). 
The uncontrolled motion is represented by an equation of 
the fifth degree, of which ofie r of the roots is zero, signi- 
fying. that the airplane is insensitive to direction in az- 
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imuth. In the remaining quartic, one of the roots is very 
large compared with the remaining three, and is very near- 
ly equal to the damping in roll. -.That ia to say, the re- 
sistance of the airplane to rate of roll is. so largo that 
the rolling ' subsidence preserves its character in spite of 
the othor lateral motions of the airplano. The rolling 
motion disappears almost at once in normal flight and may 
give trouble only 'at the stall. 

Of—the remaining throe, roots, two are conjugate com- 
plex and one is real, defining an oscillation and a subsid- 
ence (or diverg-ence ) . The only airplane "characteristics 
influencing the roots at the disposal of the designer are 
the amount of dihedral and the amount of fin and rudder 
area. Generally speaking, a large amount of static sta- 
bility causes the real root to become negative (spiral in- 
stability) while an excoss of negative static stability 
beyond a small minimum value will cause the oscillation to 
become first of increasing magnitude and then to separate 
into a rapidly increasing exponential mode (reference 36). 



Controlled Lateral Motion 

The -lateral motion is controlled by the rudder about 
the. yawing axis and by the ailerons about the rolling axis. 
The ailerons, as a rule, in addition to exerting a rolling 
moment, will also apply a yawing moment, usually of oppo- 
site sign to the direction of the desired turn and of the 
order of a tenth the rolling moment. 

The rudder .control , if actuated by azimuth indica- 
tions, makes up for the natural deficiency of the airplane 
in azimuth. But the provision of a sense of direction 
does not guarantee adequate damping of the motion. The 
use of an angle of yaw control is in this respoct not so 
fortunate as the addition of an angle of pitch control for 
the longitudinal motion. The rolling subsidence remains 
very rapid, the original short-poriod oscillation is sen- 
sibly unchanged ( reference ' 6 ) ,'• but the slow spiral diver- 
gence or subsidence formerly present has now become a long- 
period banking and yawing oscillation which may be poorly 
damped and which depends on the dihedral of the- airplane 
for the regulation of its period, of the order of 15 to 20 
seconds (reference 27). 

The motion Can be modified by operating the rudder or 
the ailerons according to other characteristics of the mo- 
tion. 
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Law of Operating Control 

Hans (reference 37) has given the following table of 
disturbance detectors, according to the indications of 
which the ailerons or rudder can "be operated. 



1. Vane with vertical 

axis 

2. Free gyroscope 



3. Free gyroscope 



4. G-yro scope produc- 

ing precesaional 
■couple 

5. Gyroscope, or dif- 

ference in lin- 
ear speed of wing 
tips 

6 . Pondulum in ZOY 

plane , or acoeler- 
omotor along OY 

7,8 Torsional accoler- 
omoters about 
X and Z axes 

9. Compass 



Recorded quantity 
Angle of sideslip 



Yaw with respect to 
axes fixed in 
space • 

Roll with respect to 
axes fixed in 
space 

Angular velocity of 
rolling 



Angular volocity of 
yawing 



Direction of appar- 
ent gravity 



Angular accelera- 
tion about OX and 
OZ 

Yaw with respect to 
earth's magnetic 
field 



Symbol 
v/ u 

CP 



g sin op 
+dv/dt+Vr 



It may be stated as fundamental that the primary pur- 
pose of tne latoral control is to give tne airplane sensi- 
tivity in azimuth. In this respect, it differs from the 
longitudinal control that operated to improve course sta- 
bility already present. Secondary control components are 
then used to improve the resulting motion. 
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The Sperry control operates the ailerons according to 
angle of bank to improve the motion, and the Smith control 
employs both aileron and a component of rudder motion pro- 
portional to angle of- Dank. The Askania control mores the 
rudder according to "both angle and rate of yaw. Garner 
has shown (reference 6) that all these secondary controls 
affect the damping of the oscillations. The rate of yaw 
control, hoivever , affects the short oscillation almost ex- 
clusively instead of the long-period "course" oscillation 
which it was intended to improve. The aileron controls 
improve the motion "by modifying the distribution of damp- 
ing although they are simple displacement controls. 



Lag in the Control 

There are two separate controls, as a rule, for the 
rudder and the ailerons. Consequently, two additional 
equations of-motion are introduced, one for each new do-, 
gree of freedom, and the determination of the complete mo- 
tion will require the solution of a ninth-order equation. 

It is possible that the highest order terms may "be 
neglected, or that lag in the aileron control may bo nog- 
lectod compared with lag in the rudder control, because of 
the much greater damping in roll than in yaw. Sufficient 
calculations are not yet available, however, to endorse 
such simplifications. 

Following the method of introducing control inertia 
and damping that was treated in deta.il under longitudinal 
control, the complete equations of controlled lateral mo- 
tion will now be presented, for an airplane with rudder 
and ailerons moved according to angles and first and sec- 
ond derivatives of yaw and roll. This arrangement is 
analogous to the control provided by the Sperry Gyropilot , 
except for the addition of the derivative components. 
Controls using. other laws of operation can be similarly 
analyzed. 

Equations of Motion 

For simplicity, the airplane will be considered in 
level flight -so that 8 q = o , and the product cf inertia 

E will be considered small enough to be neglected. The 
equations of lateral motion are thon 
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Y = m(dv/dt+Ur) (70) 
L = A dp/at (71) 

¥ 

H = 0 dr/dt (72) 

Where Y is lateral force, L is rolling moment, and 
N is yawing moment. 

There are two separate controls, one for the ailerons, 
and one for the rudder. Let 

3? , force on aileron control. 

T , force on rudder control. 

t, displacement of aileron control. 

Tt, displacement of riiddor control. 

Force derivatives for the controls are ohtainod in the 
same manner as force and momont derivatives for the air- 
plane . . Thus , 



F | = (^control) OF/afp "to. „ . (73) 



Then the free motion of oach control, with no forc- 
ing function or extornal force acting can "be written as 



Aileron control (D + DFg + Eg ) I = 0 (74) 

Rudder control (D 2 + DTfj + i\ = 0 (75) 

These two equations define the natural period and 
the damping of each control. 

The full set of simultaneous equations for the com- 
plementary solution of the motion of the controlled air- 
plane is given hy — 

Side force 

(D-Y v ) v + (-Yjp) cp + (DU) + 0 + 0 = 0 (76) 
Rolling moment 

(-L T v + (lf-DL p )q+ (~DL r )^ + (-L|)|0 ' (77) 
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Yawing moment 

(-N v )v + (-DNp)cp + (D 2 ~DN r )^ + (-Ng )| + (-N^)^ = 0 (78) 

Aileron control force 

0 + (-D 2 i^-JDPp-i^ )cp + 0 + (D 2 +DF^+P^)t +0 = 0. (79) 
Rudder control force 

0 + 0 + (~D a Tty-.DT r ~Oty W + 0 + (D a +DG?j]+(Pf| )f| = 0 (80) 

Dimensions of Determinant 

It is again desired to have the derivatives dimon- 
sionless. The dimensions of" the determinant resulting 
from the foregoing equations are 



T - 1 


LT _S 


LT~ S 


0 ■ 


0 




T~ 3 


T~ 3 


IT" 1 T - 3 


0 


t-1-,-3 


■ T~ 3 




If" 1 T _ 3 


L -i T -a 


0 


LT~ 3 


0 


T' 3 


0 


0 


0 


LT~ 8 


0 ' 





In order to make this expression dimensionless , mul- 
tiply 

the first line by T 

the second and third lines "by TL 

the second and third columns by T/L 

the fourth row and column by T 

the fifth row and column by T 

The unit of time is again 3? ~ m/(p/2 SU), but the 
unit of length according to Koppen.'s notation is the wing 
semispan (b/2). 
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Expansion, of the Determinant 



As in the longitudinal case, the determinant can "be 
expanded in terms of ideal control derivatives and the 
control frequencies and the damping ratios, so that it 
will not he necessary here to evaluate the coupling terms 
L | » ^| i K-pj, etc . 

Certain minors of the full determinant will "be used 
sufficiently often to justify a general symhol. Let the 
dimensionless determinant for the uncontrolled airplane 
be A.. Then 



4 o =s 



Also let 



d - 



- n. 



d a - dl. 



-dn, 



dp. 
— dt r 
d a - da. 



(82) 



d - y v 
-I* 



dp 



(83) 



d - y, 



-a. 



dp 
d a - dn. 



(84) 



A 3 = 



d - y T 



d a — dl, 



Now also assume the shorthand notation that 



(85) 



....■jJa + B1 p» + 1 'p 



(86) 



D n$ + Dn pc + n^ 



^ D a (a? na /2Tra?) a + D2£ 2 (T na /2TTT) + 1 



a 1 na< 



(87) 
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D s n& + Dn ro + jty 

n TIT = (88) 

D a (T ni /27T0?) a + 1)2^(^^2711) + 1 

Note that, if the yawing moment of the ailerons is 
assumed to "be proportional to the rolling moment over tho 
normal range of aileron angles, 

ngTj. = k .Igy (89) 

where k is a constant of proportionality. 

In terms of these minors and fractions, the fifth- 
order determinant can now bo expanded into a form more 
convenient for treatment. 



Expanded Form of the Determinant 

The expansion is of the form 

uncontrolled aileron yaw aileron roll u rudder yaw 
airplane com- component component + component 

ponent 

rudder yaw 
, aileron roll _ 

+ side force = 0 " (90) 

component 

Written symbolically , . this equation becomes 
A o + n f F A i - Igp A a - n T]T A;5 + n T}3? . £J- ( D ~y v ) =0 (91) 

It is interesting to not© that, whon the rudder is 
moved according to tho angle of roll, the effoct corre- 
sponds to the yawing action of the ailerons and is there- 
fore included in the second term of equation (91 ), Thus, 
rudder movement according to this relation can be made to 
balance exactly tho adverse yawing moment of the ailerons. 
Tho operation is ono that human pilots perform instinctivo- 
ly. 
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Effect of Ideal Controls 



It is interesting to compare the effects of the var- 
ious components of control outlined in the foregoing on 
the motion of the airplane. If the simple displacement 
controls of the rudder and ailerons are perfectly quick 
and powerful, the airplane is completely restrained in roll 
and yaw, and equation (91 ) reduces to simply 



D - y v = 0 



(92) 



which denotes a subsidence in sideslip, determined by the 
side force produced by sideslipping. 



minant 
form 



The uncontrolled motion is represented "by the deter- 



A Q which, when expanded, gives an equation of the 



D 5 + bE 4 + cD 3 + dK 8 + eD = 0 



(93) 



The effect of the control derivatives is to increase, or 
decrease these coefficients or to raise them to higher pow- 
ers of D. Suppose, first, that the controls are ideal, 
and that they are simple displacement controls with no de- 
rivative components. Then the additions to the coeffi- 
cients can he tabulated. 



Add to 
coeffi- 
cient 



Rudder mo- 
ment pro- 
portional 
to yaw 



-Hn llf (L p y v ) 
-iuLty.(iil T 0 L '). 



Rudder moment 
proportional 
to angle of " 
bank 



0 

0 .. . 

-Hn^dj. ) 
IWq, (l r .y v +y.i v > 
0 



Aileron yaw- 
ing moment 
proportional 
to angle of 
bank 

0 
0 

-Hn^Uj, ) 

*"Vd r y v +nl v ) 
o 



Aileron rolling 
moment propor- 
tional to an- 
gle of bank 



0 

p,l C p(n r +y v ) 
-Hl (?) (n r y T + ^n r ) 

. 0 



(94) 



When "both rudder displacement' proportional to yaw and ai- 
leron displacement proportional to angle of bank are pres- 
ent, there is an additional cross-product term" 
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Hn^ nlq,(D - y v ) (95) 
Discussion 

In. order that there "be course stability, f must he 
greater than zero. If a rudder control is used alone, 
then for f to he positive l T must ho greater than zoro 
numorlcally (negative in this system of units; l v is the 

derivative depending on dihedral and with l v zero (slight 
negative dihedral.), the pilot .cannot maintain course sta- 
bility without the aid of the ailerons. Addition of tho 
simple aileron control" adds a term from (95) proportional 
to y v , the side force in sideslip, and this term in- 
creases the course static stability. 

If f were z ero, e would determine the spiral sta- 
bility of the airplane. For most airplanes, e is nega- 
tive, indicating a divergence. If a simple aileron con- 
trol is added alone, the positiveness of e can be in- 
creased and the airplane made spirally stable regardless 
of. a small deficiency of dihedral. 

When the rudder control is added, the airplane cannot 
be spirally unstable, and the term e becomes indicative 
of, but not the criterion of, the damping of the course 
oscillation. The damping is evidently improved by large 
dihedral., side area, and control derivative l^. 

None of those simple controls affects b, which very 
closely ropros,onts the rolling subsidence. Therefore, 
this component of the motion will be substantially un- 
changed by the controls in normal flight. 

A satisfactory approximate factorization has not yet 
been developed for the short oscillation, but its damping 
is probably improved by a large value of - coefficient c. 

The effect of the first and. the second derivative 
components in the controls is to raise the individual 
terms of (94) by one and two rows, respectively, and to 
add them to the original simple displacement components. 
That is, in tho first column, the n^ control contributes 
-M.n\|; to the c coofficiont. An n rc component will con- 
tribute -(j,n rc to the b torm.' And so on with tho re- 
maining terms. 
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Moro detailed conclusions will require extensive cal- 
culations. • - 

Nonideal Controls 

There is no difficulty in expanding equation (91), 
when the controls do not have negligible inertia. Neglect- 
ing the higher derivative components of the control, which 
are treated in a similar manner, and letting 

D 2 (T ai /2TrT) 3 + D2 5 1 (!P ni /2TTa?) = ■ C x , etc. <96) 
Equation (91 ) expands to 

(Oi + DCCg+l) A Q + (C x +l) nqjAi - (Oj + l) l.yA s - 

- (C 2 +l) nv|/A 3 + n^cp^-SV ) = 0 (97) 
or, when separated into ideal and nonideal components:' 
Ideal A Q + n^Aj - l^A a - n^Ag + n^p-y.^) 

First order contribution + ( C 1 +C 2 ) A 0 +Ci (n<p ^x^cp As )-C a n^A 3 
Second order contribution + C x C 3 A Q => 0 (98) 

T7hen either the rudder or the aileron control is being 
considered alone, the frequency expression simplifies con- 
siderably. In its full form as given, it involves a ninth- 
order equation, but Of Ga is probably very small compared 
with the first-order terms, and possibly may be neglected. 

The stability of the motion can then be investigated 
by means of the natural-frequency expressions for the sep- 
tic. 



January 18 , 1939 . 
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Figs. 2,3 




Figure 3,- Effect of displacement control on period and damping 
of phueoid. 



